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Abstract. Scalar cosmological perturbations are investigated in the framework of a model 
with interacting dark energy and dark matter. In addition to these constituents, the inho¬ 
mogeneous Universe is supposed to be filled with the standard noninteracting constituents 
corresponding to the conventional ACDM model. The interaction term is chosen in the form 
of a linear combination of dark sector energy densities with evolving coefficients. The meth¬ 
ods of discrete cosmology are applied, and strong theoretical constraints on the parameters 
of the model are derived. A brief comparison with observational data is performed. 
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1 Introduction 

The coincidence problem is among the greatest challenges in modern cosmology. It consists in 
the observational evidence that within the conventional ACDM model the energy fractions of 
dark energy and dark matter are comparable (or, in other words, represent quantities of the 
same order) in the Universe today, despite the fact that these constituents evolve in a quite 
different way; namely, the A-term does not evolve at all, while the energy density of cold dark 
matter (CDM) decreases during the cosmological expansion. Dark energy represented by the 
A-term completely dominates the evolution of the future Universe. However, at the early 
evolution stage it is the CDM that plays the leading part. In this connection the following 
natural question arises: why are the dark sector contributions comparable today? Rejecting 
randomness as an explanation for such a coincidence, one can alleviate this challenge by 
introducing direct non-gravitational dark energy - dark matter interaction [1-5] (see also the 
important early papers [6, 7] on coupled scalar fields and the detailed review [8] on dark 
energy dynamical behavior). From the particle physics point of view, the case of uncoupled 
dark energy and dark matter is anyway a particular one, and in general their coupling may be 
nonzero and natural [5] . As an example of such a coupling, let us mention the mass of CDM 
particles depending on the dark energy field [9] . It should be mentioned in addition that the 
introduction of dark sector interactions is also reasonable for the holographic scenario [10, 11]. 
Besides, the data of observations from the Abell Cluster A586 [12-14], measurements of the 
baryon acoustic oscillations of the Lyman alpha forest from high redshift quasars [15] as well 
as the cosmic microwave background and redshift-space distortions measurements [16] may 
indicate the reality of these interactions. 

In the present paper, we use the methods of discrete cosmology devised in [17-19] 
for the ACDM scenario to investigate a cosmological model with dark energy - dark matter 
interaction. Our approach consists in describing scalar perturbations at the stages of the cos¬ 
mic evolution when inhomogeneities such as galaxies and galaxy clusters have already been 
formed by representing the nonrelativistic pressureless matter (dust) as a system of separate 
gravitating point-like particles. This is in accordance with the observational evidence that the 
typical distances between separate cosmic bodies essentially exceed their dimensions. This 
well-grounded idea turns out to be interrelated with gravitational potentials and Newtonian 
equations of motion on the FLRW background, concepts which are commonly used in mod¬ 
ern N-body simulations. At the same time, it leads to theoretical as well as experimental 
restrictions for numerous extensions of the conventional model and its alternatives. For ex¬ 
ample, methods of discrete cosmology have already been applied successfully to cosmological 
models including a perfect fluid with a constant negative parameter in the linear equation 
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of state (e.g., quintessence, the phantom field, topological defects such as cosmic strings and 
domain walls) [20], Chevallier-Polarski-Linder (CPL) and other linear parametrizations of 
the equation of state parameter [21], quark gluon plasma that has escaped hadronization and 
survived up to now in the form of quark-gluon “nuggets” [22, 23], nonlinear gravity [24], and 
the Lattice Universe topology [25]. 

Here, we use the same approach with respect to the inhomogeneous Universe containing 
the standard ACDM constituents as well as interacting dark matter and dark energy. We 
start by analyzing the equations of the scalar perturbations theory for the nonzero interaction 
term, then derive the conditions which must be satisfied by its parameters, and conclude by 
briefly comparing the established theoretical limitations with the observations. 


2 Dark sector interactions and scalar cosmological perturbations 


First of all, let us write down the standard Friedmann equations describing the evolution 
of the homogeneous background in the framework of the conventional ACDM model supple¬ 
mented with the dark sector interactions: 


3 


kTq -|- a -|- K££)M + + E,£rad j 


( 2 . 1 ) 


2^' + n‘^ + K, 


= A - KpDE - E,Prad ■ 


( 2 . 2 ) 


Here, T-L = a'/a; the prime denotes the derivative with respect to the conformal time 77 ; 
fC = 0, ±1 represents the spatial curvature; «(??) is the scale factor; k = 8ttGn/c^, where Gjy 
is the Newtonian gravitational constant, and c is the speed of light; Tq = p(?/a^ is the average 
energy density of the noninteracting nonrelativistic matter (both dark and visible), where p 
is the corresponding average comoving rest mass density; A represents the cosmological con¬ 
stant, while e_D_E(7?) are the average energy densities of the (generally speaking) 

spacetime-dependent dark matter and dark energy, respectively, interacting with each other; 
Pde is the pressure of the interacting dark energy (the pressure poM of the interacting dark 
matter is assumed to vanish, since it is also assumed to be completely nonrelativistic, similar 
to the noninteracting one); finally, £rad{'n) is the average energy density of radiation, and 
Prad = ^radl‘^> is the Corresponding radiation pressure. 

For simplicity and illustration, we focus our attention here on the vacuum-like equation 
of state pde = —£de for the interacting dark energy, similar to the noninteracting one 
represented by the cosmological constant A (see also, e.g., [26, 27], where the same choice is 
made). In this case, one can also use the notation k£de = A, and we can write npoE = ~A. 
Here, evidently, h.{ri) is not a constant; it varies with time in view of the presumptive non- 
gravitational interaction. 

In the framework of the discrete cosmology (mechanical) approach to cosmological prob¬ 
lems inside the cell of uniformity, developed recently in [17-19] (see also the recent papers 
[28-32] on related issues), let us consider the following FLRW metric perturbed by inhomo¬ 
geneities: 


ds^ 


(I -|- 2^)d'rf‘ — (1 — 2^)jai3dx'^dx^ 


a, /3 = I, 2,3 . 


(2.3) 


Here, ^{r],r) = p{r)/ (c^a), where the introduced function p depends only on comoving 
spatial coordinates r and does not depend on time within the adopted accuracy (both the 
nonrelativistic and weak field limits are applied), in agreement with [33-35]. Really, this is 
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a general solution of the equation <!>' + = 0, where contributions from peculiar velocities 

of inhomogeneities are neglected. A noteworthy feature of this solution lies in the fact that 
for a single delta-shaped inhomogeneity it has the correct Newtonian limit $ ~ 1/i? when 
~ 1/r, where R = ar \s the “physical” distance (see also [17, 18] for details). The scalar 
perturbation $ satisfies the following system of linearized Einstein equations [36, 37]: 

A<h -|- 3/C<h = 2^®^ {^'^0 T + ^^rad) j (2-4) 


+ + 'H^) {6pdE + ^Prad) = {^^£DE + ^SSrad^ , (2.5) 


where A stands for the Laplace operator defined with respect to the spatial metric coefficients 
7 q,^, and the fluctuation of the energy density Tg of the noninteracting nonrelativistic matter 

0 _ 6p{r)c^ 3pc^$ _ <5p(r)c^ 3p(/?(r) 

01 Q ' A 

a'^ a'^ 

Here, (5/3(r) stands for the fluctuation of the rest mass density in comoving coordinates. From 
(2.5) with the help of (2.1) and (2.2) we immediately obtain 


( 2 . 6 ) 


^ + eoMir])^ + ^Sradiv)^ = SedE - bsrad • (2.7) 

c^a 3 Ao 3 

It is worth noting that this equation agrees with the corresponding equation (40) in [38] 
in the absence of interacting dark constituents. The third term on the left-hand side (Ihs) 
of (2.7) behaves as 1/a®, which is outside the limits of our adopted accuracy and therefore 
should be dropped, while the first term on the Ihs as well as the second term on the right 
hand side (rhs) both behave as 1/a^. In other words, the background radiation contribution 
can be disregarded. Then we get from (2.7): 

5eDE = ^ + EDMiv)^ + berad ■ ( 2 - 8 ) 

c^a 3 

The dependence of the second term on the rhs of (2.8) on the scale factor a is different 
from l/a"^, since the quantity edm{ii) does not behave as 1/a® (this behavior is inconsistent 
with the dark sector interactions). Moreover, this term cannot be neglected because the 
cosmological model under consideration includes interacting dark components with the energy 
densities comparable with the energy densities of noninteracting constituents at present, as 
well as in the (at least) near past and in the (at least) near future. Consequently, even if 
^Erad = ~3p99(r)/a'^, and the hrst and third terms on the rhs thus exactly compensate each 
other, we would still have 6ede 7^ 0; the interacting dark energy thus cannot be homogeneous, 
in contrast to the noninteracting one represented by A = const. 

The substitution of (2.8) into (2.4) leads to the following equation: 


A(/?(r) -F 3)C(f{r) = (^6p{r)c^ + + u^Sedm + a^EDMiv)^^^ + b^^Erabj , 

(2.9) 

whence 

/Sip{Y)+3K,p{r) = [<5/0(r)c^ + /(r)] , ^EdM = + • 

Zi Qj C CL O CL 

( 2 . 10 ) 
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Here, the introduced function / depends only on comoving spatial coordinates r similar to 
the function (p, while its dependence on time should be disregarded. 

It is worth mentioning that, as one can easily verify, the derived fluctuations (2.8) and 
( 2 . 10 ) satisfy the energy conservation equation 


+ ^T~i{5£tot + Sptot) + i^tot + — — 0 , 


( 2 . 11 ) 


where the total energy density and pressure of the interacting dark components read, respec¬ 
tively, etotiv) = £de + £dm and ptotiv) = Pde + Pdm = -^de- 

The supposed interaction between the dark constituents of the Universe is usually de¬ 
scribed by the following equations [1-5, 9, 12, 13, 39-42]: 


^'de T ^T^{£de + Pde) — ^'de ~ ~Q > 

^'dm + ^'H{£dm Pdm) = ^'dm + ^T~L£dm = Q ■. 


( 2 . 12 ) 

(2.13) 


where the phenomenological interaction term Q often represents a linear combination oi £de 
and £dm with some (generally speaking) time-varying coefficients Ti and r 2 , which depend 
neither on £de, nor on £dm, and may change in time only because of some postulated 
reaction rate evolution of unknown nature: 

Q = ^i£de+ ^2£dm ■ (2-14) 

Consequently, the fluctuations must satisfy the appropriate perturbed equations: 

^^'de T 3'H((5e£)E -|- 5pde) + {^de + Pde){~^^') = ^^'de ~ ~^Q ! (2-15) 

^^'dm + ^T~l-{d£DM + ^Pdm) + if DM +Pdm)(~3$^) = + ZT-LdiDM + £dm(~3<I>^) = 5Q , 

(2.16) 

where 

5Q = Ti6£de+ ^25£dm ■ (2-17) 

Hereinafter, we focus our attention on the most popular case of the following parametriza- 
tion [1-5, 9, 13, 39-42]: 

ri=777i, r2 = 7772, (2-18) 

where 71 and 72 are constant (time-independent) parameters. It is important to stress 
that the quantity % is treated here as characterizing the global average expansion rate. 
Consequently, neither % nor ri ^2 are perturbed when passing from (2.14) to (2.17). At the 
same time, according to the reasoning presented in [43] (see also [44, 45]), one should not 
ignore the perturbation 5H describing the difference between the local expansion rate and 
the global one, which is defined by means of the relationship 


< « -{n + 5 %), 

’ a 


(2.19) 


where i = 0,1,2,3, stand for the 4-velocity components. However, it is easy to demon¬ 
strate that within the adopted accuracy dPL = 0. Really, since we neglect the peculiar motion 
of inhomogeneities, the only nonzero component of the 4-velocity is ~ (1 — 4>)/a. Then 
the calculation of the covariant derivative on the Ihs of Eq. (2.19) on the basis of the metric 
(2.3) gives 


u.. 


1 - 4 > 


-7 


4a' 1 - 4> 24>' 


m _ 377 4 > 

a 


a 


34>' 

a 


( 2 . 20 ) 
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up to the first order of smallness. Further, taking into account that <h(r/, r) = (/j(r)/ (c^o), 
we conclude that ~ 3'H/a. Then the comparison with (2.19) immediately gives 51-1 = 0. 
Thus, disregard of the perturbation dH is really justihed within the accuracy of our approach. 
Substituting (2.17) and (2.18) into (2.15), we obtain 

^^'de = + 125£dm) ■ (2-21) 


Furthermore, substituting (2.8) and (2.10) into (2.21), we get after some computation 

^2 


— 2 2 2, -p f \ 

lieDEii) = (4-7i + 472)^ + (4-7i)ggAf(?7)+ (4-7i+472)ferad-72 A • 

a-^ 3(/?(r) a^(^(r) 


In view of (2.11), Eq. (2.16) gives nothing new and may be dropped. The substitution 
of (2.14) and (2.18) into (2.12) and (2.13) leads to the equations 

=— {ii£de + 12 £dm)-, = -[jiSde + (72 ~ 3)g_DM] ■ (2.23) 

da a da a 

After some computation with the help of (2.22) and (2.23), one can show that 

— 2 2 2 p / \ 

- (4-71+ 472)gDM(h) = (4 - 71 + 472)^ (4 - 7i + 472 )fer-ad " 272 . . ( 2 . 24 ) 


3ip{r) 


‘ a‘^ip{r) 


First, let us consider the case 4 — 71 + 472 = 0, which gives 72 /(r) = 0. If 72 = 0, and 
thus 7 i = 4, there is a contradiction with Eq. (2.22). Therefore, 72 7 ^ 0, while /(r) = 0. 
According to (2.10), in this case the nonrelativistic gravitational potential is determined only 
by inhomogeneities of noninteracting matter as in [17-19]. From (2.22) and (2.23) we obtain 


£DE{'n) = Aa 3 ( 1 + 72 )^ £DM{ri) 


^1 + —^ Aa 3 ( 1 + 72 )^ ^ _ const. 


(2.25) 


If we demand in addition that both energy densities in (2.25) are positive, then A > 0 
and —1 < 72 < 0, so £de{v) EDMiv) decrease in the expanding Universe in view of the 
inequality 3(1 + 72 ) > 0. According to (2.8) and (2.10), the same statement holds true for 
their fluctuations. Besides, Q > 0, and this means an energy transfer from dark energy to 
dark matter. Exactly this transfer direction is required by the second law of thermodynamics 
as well as for solving the coincidence problem [3, 13, 40, 46-48]. It is worth noting that if 
the parameters 71 and 72 are supposed to be equal, then we have 71 = 72 = —4/3, so the 
condition 72 > —1 does not hold true. 

In the second case we have 4 — 71 + 472 7 ^ 0. Therefore, from (2.24) it immediately 
follows that 


^£-rad 


3w(r) 


£DM{r]) = 


272 cV(r) 

4 - 7 i + 472 a‘^(p{r) ' 


(2.26) 


Here, we take into account the above-mentioned fact that the quantity £dm{ii) does not 
behave as l/o3, so the only possibility to compensate the last term on the rhs of Eq. (2.24) 
is to require that £dm{ii) ~ 1/®^ and 72 /(r) 7 ^ 0. From (2.22), (2.23) and (2.26) we get the 
equality 2 — 71 — 272 = 0 and the expressions 


SDMiv) 


72 c ^/( r ) 

1 - I - 372 a'^(p{r) ’ 


- ^ ^ 72 cPf{r) 

^BEin) = jEDM = 2(1 + 3b,) ■ 


(2.27) 
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As in the previous case, both energy densities eDEiv) eDM(^) decrease in the 
expanding Universe together with their fluctuations. Besides, we again demand additionally 
that £de,dm > 0 and come to the conclusion that the energy transfer from dark energy 
to dark matter takes place, as before, in view of the inequality Q > 0. It is also worth 
mentioning that if the parameters 71 and 72 are equal to each other, then their common 
value is 2/3. 

Expressing the function /(r) from (2.27) and substituting the result into (2.10), we 
derive the following equation for the nonrelativistic gravitational potential; 


A(/?(r) + 3IC(p{r) - X(p{r) = ^Kc^6p{r) = 47 rGAr[/ 9 (r) - p], A = ^'t ko^edm( il) = const. 

2 272 

(2.28) 

For simplicity, let us restrict ourselves to the case of flat spatial geometry (/C = 0) and 
require that A > 0, so either 72 < —1/3 or 72 > 0. Then the solution of Eq. (2.28) for a 
system of gravitating masses m* with the comoving radius-vectors rj, characterized by the 
rest mass density 

p{r) = ^ mi6{r - r*), (2.29) 


reads 


<^(r) = -GnYI 


rrii 


r - r,; 


exp 


-y/X\r - Til) -h 


AttGnP 

A 


(2.30) 


It is very important that this function is characterized by the zero average value (p = 0 
for any mass distribution, in contrast to the potential determined within the bounds of the 
conventional ACDM model, which may have a nonzero value after performing a procedure 
of spatial averaging (for an example of the corresponding mass distribution, see [17, 18, 49]). 


3 Comparison with experimental restrictions, and generalization 

Summarizing, we have arrived at the following two allowed cases: 

I. 4-71 -1-472 = 0, / = 0, —1<72<0, A>0 

II. 2 - 71 - 272 = 0, /72(1 372 ) > 0 

Provided that the most widespread case of flat spatial geometry (/C = 0) is studied, the 
second class becomes even more constrained: 

II. 2 - 7 i - 272 = 0, / > 0, 72 < -1/3 or 72 > 0 

Now let us illustrate these constraints with three popular examples: 

A. Q = 'HjiSde (72 = 0) 

B. Q = 'H'y 2 £DM (71 = 0) 

C. Q = 7^7 {sde + edm) (71 = 72 = 7 ) 

One can easily verify that the model A satisfies neither the first set of conditions (I), 
nor the second one (II). The models B and C also do not satisfy the conditions I, while the 
conditions II lead to the following allowed interaction terms, respectively; Q = Hedm and 
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Q = {2/2>)'H {ede + £dm)- If the noninteracting dark sector components were absent, both 
expressions would be incompatible with known observational restrictions for the parameters 
7 i and 72 , see [9, 13, 39-41]. Consequently, two of the above-mentioned interaction terms are 
consistent with the theory of scalar cosmological perturbations and with observational data 
only if in addition to interacting dark constituents of the inhomogeneous Universe there are 
noninteracting constituents weakening the empirical limitations. Actually, this is the reason 
why we have included the noninteracting dark sector components (in particular, the constant 
A-term) in the model. In other words, after imposing theoretical constraints, the interacting 
dark energy does not necessarily ensure, e.g., the late-time acceleration of the Universe 
expansion (in this case the “dark energy” notion itself loses the initial sense). Therefore, 
the constant A-term may be still needed for reaching agreement with the supernovae and 
other data. Thus, it is reasonable not to exclude such possible contributions from the very 
beginning. 

The finding that three popular interaction terms quoted above are in general incompat¬ 
ible with the derived constraints actually originates from Eq. (2.22) and is certainly valid for 
the investigated simplest case of the vacuum-like equation of state pde = —£de- Tbe more 
general and complicated case of an arbitrary constant parameter u in the linear equation of 
state Pde = ^£de lies beyond the scope of our current investigation, representing a subject 
of a separate forthcoming paper. However, let us laconically generalize Eq. (2.22) for uj < —1 
or —1 < w < 0. First, instead of (2.8) and (2.10) now we have, respectively, 

Sede = - SDMiv)-^^ - “ (1 + w)ez)£( t/), (3.1) 

ioa^ ioc^a 3io ioc^a 

- + {1 + uj)ede{v)^^^ ■ 

ujc^a 
(3.2) 

Second, the background energy densities £DE{fl) ^DMip) now satisfy the equations 

= —[(71 + 3-1- 3lo)£de + 12 £dm], = -[jl^DE + (72 “ 3)£dm] ■ (3.3) 

da a da a 

Finally, the generalization of Eq. (2.22) reads: 

[71 - (1 + a ;)72 - (1 + w)(l -F 3uj)]£DEi'n) 

- 2 

= [1 - 3a; - 71 -h (1 - 3 a;) 72 ]^ -F [1 - 3a; - 71 -F (1 + u})'y2)]£DM{r]) 

+ [1 - 3a; - 71 -h (1 - 3 a;) 72 ] ^ ^. dErad + ^ 72 % 4 ^ • (3.4) 

3(/9(r) a^(p{r) 

The derived condition (3.4) should be taken into account along with Eqs. (3.3) and 
plays the role of the additional restriction for the set of three initially free parameters uj, 71 
and 72 . 


Sedm = — 1 3 — 


UJ 


/. 1 f 


UJ 


c^a 


-H3 + 

UJ 


^£rad d 


/(r 


4 Conclusion 

In the present paper, we have investigated the scalar perturbations in the framework of the 
cosmological model with dark energy - dark matter interaction. The inhomogeneous Uni¬ 
verse is supposed to be filled with the standard noninteracting constituents of the conventional 
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ACDM model as well as with interacting dark matter and dark energy (with dust-like and 
vacuum-like equations of state, respectively). The interaction term is taken in the frequently 
used form (2.14) with the coefficients (2.18). We have applied methods of the discrete cosmol¬ 
ogy approach in order to impose theoretical constraints on the model under consideration. 
These constraints originate from the derived condition (2.22) and are generally incompatible 
with the experimental data. Then we have generalized this condition to the case of an ar¬ 
bitrary parameter u in the linear equation of state of interacting dark energy. As a result, 
it takes the form (3.4) and can be used again for restricting the parameters entering the 
interaction term, as well as the equation of state parameter u itself. 
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